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DYNAMIC  RESPONSE  OF  A  CONSTRAINED  FIBROUS  SYSTEM 
SU8JI  'TED  TO  TRANSVERSE  IMPACT 

PART  II:  A  MECHANICAL  MODEL 
ABSTRACT 

A  "finite  element"  mechanical  model,  capable  of  utilizing  experimentally 
determined  material  properties,  is  presented  for  analyzing  the  deformation 
processes  in  a  finite,  constrained  flexible  filament  subjected  to  transverse 
impact.  ” 

The  model  is  then  used  to  analyze  the  results  of  experimentally  (trans¬ 
versely)  impacted  finite  lengths  of  nylon  yarns  rigidly  clamped  at  both 
ends.  The  experimental  analysis  is  presented  in  Part  I  of  this  work. 

The  analysis  is  used  to  determine  the  effects  of  an  elastic  foundation 
transversely  supporting  a  flexible  filament  subjected  to  transverse  impact. 

Finally,  results  of  a  finite  element  analysis  of  a  transversely  impacted 
isotropic  membrane  are  presented  to  indicate  the  direction  for  future  re¬ 
search. 

An  appendix  is  attached  which  demonstrates  the  ability  to  include  rate- 
dependent  material  properties  in  the  analyses  described  above. 
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INTRODUCTION 


The  transient  deformation  of  constrained  fibrous  systems  subjected  to 
high-speed  transverse  impact  is  being  studied  to  gain  a  better  understanding 
of  the  performance  of  textile  and  felt  materials  used  for  armor  systems.  A 
description  of  the  experimental  portion  of  this  program  is  given  in  Part  1 
of  this  investigation.1  Although  there  has  been  fairly  extensive  research 
into  the  properties  of  individual,  unconstrained  fiber  components  with  trans¬ 
verse  impact  at  ballistic  velocities,2"5  there  has  been  a  singular  deficiency 
in  defining  how  these  properties  translate  into  the  system  behavior,  i.e., 
the  constrained  fiber  including  the  boundary  interactions.  The  missing  in¬ 
gredient  has  oeen  a  mechanical  model  which  can  utilize  the  fiber  properties 
obtained  from  a  tractable  experiment  to  predict  the  behavior  of  the  system 
under  practical  conditions.  This  is  essential  both  for  the  evaluation  of 
meaningful  armor  test  conditions  and  to  understand  the  performance  of  the 
fiber  under  in-use  conditions  as  a  component  of  an  armor  system. 

A  summary  of  the  experimental  program  in  Part  I  of  this  work1  is  given 
below.  The  "finite  element"  mechanical  model  is  then  developed  and  used  to 
predict  and  analyze  the  experimental  results  obtained  for  the  simple  con¬ 
strained  fiber  system.  More  practical  constraints  are  then  applied  to 
demonstrate  the  ability  to  perform  the  analysis  for  diverse  physical  situa¬ 
tions.  These  include  the  effects  of  an  elastic  foundation  resisting  deforma¬ 
tion  a-id  the  impact  of  an  isotropic  membrane. 

An  appendix  is  attached  which  demonstrates  the  ability  to  include  rate- 
dependent  material  properties  in  the  models  described  above. 


SUMMARY  OF  EXPERIMENTAL  RESULTS 

Tie  experimental  program  detailed  in  Part  I  of  this  work1  was  developed 
for  obtaining  fiber  properties  under  impact  conditions  and  studying  how  a 
simple  constrained  system  using  these  fibers  would  behave.  Figures  1  and  2 
illustrate  the  type  of  results  obtained  from  this  experimental  program. 

Figure  1  is  a  photograph  taken  with  several  flash  exposures  indicating 
the  gross  deformation  caused  by  a  missile  transversely  impacting  the  con¬ 
strained  system  of  a  low-tenacity  nylon  yarn  bundle  rigidly  fixed  at  both 
ends.  There  is  symmetry  about  the  point  of  impact;  therefore,  only  the  upper 
half  of  the  yarn  bundle  is  shown.  From  results  such  as  these,  for  very  short 
times  after  impact,  before  the  deformation  processes  interact  with  the  bound¬ 
aries,  a  dynamic  stress-strain  curve  can  be  obtained  for  the  fibers.4  This 
stress-strain  curve  will  be  some  average  value  over  the  rates  involved  in 
obtaining  the  data;  the  true  viscoelastic  relationships  are  considered  in 
Reference  5  and  the  appendix.  For  longer  times  after  impact,  after  at  least 
one  reflection  of  the  transverse  displacement  wave,  an  "average"  stress-strain 
curve  over  somewhat  lower  rates  may  be  obtained.3 

Figure  2  is  a  plot  of  the  energy  lost  by  the  projectile,  which  impacts  and 
breaks  the  constrained  fiber  system,  as  a  function  of  the  impact  velocity. 

The  energy  is  expressed  per  unit  mass  of  the  fiber  between  the  clamps. 
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Figure  1.  TYPICAL  PHOTOGRAPH  OF  IMPACTED  YARN  BUNDLE  SHOWING  DEFORMATION  AND  RUPTURE. 
Bundle  contain*  4  yarns  of  low  tenacity  nylon  6/6.  Impact  velocity  341  feet  per  second.  Time  of  325  microseconds 
(>e tween  flash  exposuies. 


MISSILE  IMPACT  VELOCITY  (FT/SEC) 

Figure  2.  ENERGY  STORED  IN  HIGH  TENACITY  NYLON  FIBER  SYSTEM  VERSUS  TRANSVERSE  IMPACT 
VELOCITY.'  Velocity  regime  I:  quasi-static  fiber  response  negligible  wave  propogation  effects;  material  pronerties 
from  method  in  Ref.  3.  Velocity  regime  II:  dynamic  fiber-system  response,  important  wave-boundary  interactions; 
material  properties  from  method  in  Ref.  4.  Velocity  regime  III:  dynamic  fiber  response  rate^ependeni  considerations 
become  important  and  boundary  effects  are  negligible. 
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It  is  now  necessary  to  develop  a  mechanical  model  in  order 
the  basic  fiber  properties,  which  can  be  obtained  as  mentioned 
late  into  the  fiber  system  behavior  indicated  in  Figures  1  and 


to  explain  how 

uKavo  *■  v*or»c  _ 
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THE  METHOD  OF  ANALYSIS 


Background 

The  mechanical  model,  which  appears  to  satisfy  the  requirements  for 
analyzing  constrained  fibrous  systems,  is  a  modification  of  the  "direct  analy¬ 
sis"  of  waves  in  elastic  media  developed  by  Mehta,  Koenig,  and  Davids.6’7 

The  modeling  and  subsequent  "direct  analysis"  using  physical  laws  in  an 
incremental  form,  rather  than  differential  equations,  constitute  what  has 
become  known  as  "finite  element"  analysis.  In  conventional  mathematical 
procedures,  differential  equations  are  derived  from  the  physical  laws  and 
an  attempt  is  then  made  to  solve  these  equations  subject  to  appropriate 
initial  conditions  and  boundary  conditions.  For  problems  as  complex  as 
those  considered  here,  the  differential  equations  must  be  solved  by  numeri¬ 
cal  methods,  and  even  then  the  numerical  techniques  are  often  quite  complex 
and  may  net  be  tractable. 

In  the  finite  element  analysis  the  differential  equations  are  regarded 
as  dispensable  and  the  physical  laws  are  used  directly  to  effect  changes 
in  the  system  that  adhere  to  these  laws  (e.g.,  conservation  of  energy, 
momentum- impulse  balance,  and  material  constitutive  laws) . 

The  main  advantages  of  this  analysis  are  that  it  is  easier  to  model  the 
physics  of  the  problem,  i.e.,  system  constraints,  and  experimentally  measured 
or  observable  behavior,  and  the  subsequent  numerical  procedure  is  more  di¬ 
rect  and  efficient. 

Procedure 

The  procedure  follows  the  development  of  Koenig  and  Davids.7 

The  fiber  system  is  divided  into  finite  lengths  (AL)  of  flexible  (i.e., 
no  shear  or  bending  rigidity)  filament  (see  Figure  3).  The  physical  laws 
are  then  written  in  incremental  form  for  a  typical  finite  element: 

The  impulse-momentum  balance,  relating  the  velocity  vector  (V)  and  ten¬ 
sion  per  lineal  density  (T/m)  over  an  increment  of  time  (At)  is 

Ap(j*,)=  at-  10 

The  strain  (E)  -  velocity  relationship  is 

AtU)=  tLLtlkgW  At  (z) 
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element  j 


Figure  3.  TYPICAL  FINITE  FLEXIBLE  ELEMENT 


The  tension  (]aer  lineal  density)  -  strain  constitutive  law  utilizing 
tiie  modulus  per  lineal  density  (E)  is  " 

^=E[e(J)]  -A£(i).  ( 3 ) 

Equations  (1),  (2),  and  (3)  completely  state  the  problem  of  traveling 
waves,  in  a  bounded,  elastic,  ideally  flexible  filament,  when  combined  with 
appropriate  boundary  conditions.  These  formulations  may  be  used  directly 
in  a  computer  code. 

In  order  to  insure  stability  of  the  calculations  it  is  convenient  to 
define  the  element  size-time  increment  relation  which  conserves  energy  in 
the  propagation  of  the  leading  wave  disturbance  (in  this  case,  longitudinal 
elastic  waves  propagating  at  the  velocity  C0  =  (Eq)1'2,  where  E0  is  the 
modulus  per  lineal  density  (e.g.,  grams  per  denier)  of  the  steepest  portion 
of  the  tension-strain  curve) : 

At  >  AL/Co* 

The  element  size  is  arbitrary,  yet  assures  that  a  further  reduction  in 
size  will  not  yield  any  significant  change  in  the  resulting  solution  of  the 
problem. 

The  boundary  conditions  and  physical  parameters  are  specified;  then 
equations  (1)  -  (3)  are  applied  sequentially  to  each  element  from  the 
impact  point  to  the  boundary  for  one  increment  of  time.  This  step  is  re¬ 
peated  to  yield  an  accumulated  solution  after  each  further  increment  of 
time. 

It  should  be  noted  that  no  specific  operations  are  required  tc  deal 
with  reflected  waves.  The  reflected  longitudinal  waves  are  automatically 
generated  simply  by  specifying  the  appropriate  boundary  conditions.  Simi¬ 
larly,  transverse  deflection  waves  are  automatically  genei'ated  in  order 
that  the  system  response  adheres  to  the  governing  physical  laws. 
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Elastic  Foundation 


In  order  to  simulate  the  effect  of  the  support  which  the  human  body  pro¬ 
vides  for  a  simple  fiber  system,  a  transverse  resistive  force  (R)  is  applied 
to  the  filament  proportional  to  the  transverse  displacement  (y) ,  R  =  Ky. 
Similarly,  air  drag,  R  =  D*Ay/At,  or  other  such  forces  could  easily  be  added 
without  increasing  the  complexity  of  the  analysis. 

Membrane  Analysis 

As  a  further  step  in  simulating  how  more  practical  fibrous  systems  behave, 
a  finite  element  model  for  the  transverse  impact  of  an  isotropic  membrane  (a 
first  approximation  to  a  sheet  of  fabric)  is  now  developed. 

The  membrane  is  divided  into  finite  s^mi-annuli  as  shown  in  Figure  4. 

The  physical  laws  are  then  written  in  incremental  form  for  a  typical  finite 
element.  The  treatment  of  cylindrical  elastic  waves  in  Reference  6  follows 
a  similar  development  to  the  following. 
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Figure  4.  TYPICAL  FINITE  FLEXIBLE  MEMBRANE  ELEMENT 


The  impulse-momentum  balance  relating  the  particle  velocities  (u,v)  in 
the  r,y  directions,  the  radial  (O',)  and  circumferential  stresses  and 

the  density  i />)  per  unit  volume  are  derived  by  first  considering  the  motion 
of  the  center  of  gravity  of  the  semi  annulus  ^V^CC.SO] 

in  the  x  and  >  directions,  caused  by  the  net  forces^  p  T  Fv  ^  *n  these 
directions:  *>- 

mnl()  -avtic-z  )  =  » 

m(i)  AVsic.g.)  =  ^;f<  At. 

The  velocities  of  the  center  of  gravity  must  now  be  related  to  the  element 
particle  velocities.  The  element  shoving  from  radial  position  r  to  r*Ar  has 
its  center  of  gravity  moved  from^f  to  in  time  At  thereforeAVfc(C.ST)= 

But,  the  particles  along  a  radius  of  the  element  move  from  r  to  r-*Ar 
therefore,  AU  =  ^-.  The  relation  between  the  particle  velocity  in  the  r- 
direction  and  the  center  of  gravity  motion  in  the  x-direction  becomes 
AU  =  AVVCC.gh)  .  The  particle  velocity  in  the  y-direction  is  identical 
to  the  velocity  of  the  center  of  gravity  in  the  y-direction,  AV^AtyCC.g'). 

Therefore,  the  impulse-momentum  balances,  relating  the  particle  veloci¬ 
ties  (u,v),  the  radial  and  circumferential  stresses,  and  the  density  per 
unit  volume  (where  ot  is  the  slope  of  the  membrane  with  respect  to  the  ori¬ 
ginal  plane)  are: 

a  a  ={  crru-M)  •  r(  *  vo  cosset*  i)]-oKO  r(Ocos[«(oj 
♦  05i(OCr(O-ru+i^}  2At /f>[rX<M)-rHi)}  , 

AV ={  o' (  i  *  i )  r  ( ( +J)  $in[oc(  i +1)) 

-  crr(c )  rtf)  si  n  [<x(f )]}  2  At /p[r X  i  +  J )  -  r  2(  1  >] . 


VC5) 


The  radial  (£r)  and  circumferential  )_strain-velocity  relationships 
using  the  resultant  particle  velocity  vector  V  and  radial  particle  velocity 
for  an  element  length  Afr  are: 


*  ^  V/CO-J)  -  V<i) 

A£r((;= - ^ - At  ; 


AEe«)  = 


qCQAt 

rtf)  ’ 


} 


(6) 


The  stress-strain  constitutive  lt.w  using  the  tensile  modulus  (E)  and 
Pi.isson's  ratio  (v)  is 

acrr(i)  -E—  + 

L  VAtr(i)  + 

(Computations  show  that  would  become  negative.  Therefore,  a  condition 
for  buckling  is  stipulated.  If  0$  tries  tc  go  negative,^  is  set  equal 
to  zero  in  the  constitutive  relations  above.  In  addition,  the  excess 
material  (circumferentially)  in  the  buckling  zone  (slack)  is  continually 
accounted  for,  and  reloading  must  first  take  up  the  slack.) 


Equations  (3) ,  (6),  and  (7)  completely  state  the  problem  of  traveling 
waves  in  an  elastic,  ideally  flexible  membrane  allowing  for  circumferen¬ 
tial  buckling,  when  combined  with  appropriate  initial  conditions  and  bound¬ 
ary  conditions. 

In  order  to  insure  stability  of  these  calculations  the  element  size¬ 
time  increment  relation  At  >  Aft/Cj,  was  used,  where  C0=[E<>/ /’(l  — Z''2)] 

(E0  is  the  highest  value  of  tensile  modulus  corresponding  to  the  strain 
magnitudes  encountered  ir  the  deformation  process) . 

RESULTS  OF  THE  ANALYSIS 
Transverse  Impact  of  a  Linear  Elastic  Filament 

The  following  computations  utilize  a  linear  approximation  of  the  stress- 
strain  curve  obtained  by  Smith1*  for  a  high  tenacity  nylon  yarn.  The  fiber 
studied  by  Smith  appears  to  be  quite  similar  to  the  high  tenacity  nylon  yarn 
used  in  the  experimental  phase  of  the  present  work.1 

Figure  5  illustrates  typical  computed  results  for  the  longitudinal 
strain  distribution  at  various  times  after  impact  in  a  linearly  elastic 
fiber  subjected  to  transverse  impact.  The  reflection  of  waves  at  the  bound¬ 
ary  are  automatical ty  accounted  for,  simply  by  specifying  the  fixed  boundary 
velocity  to  be  zero.  Note  that  there  is  an  interaction  oetween  the  longi¬ 
tudinal  strain  wave  and  the  transverse  filament  deflection  at  approximately 
110  psec  which  causes  a  certain  percentage  of  the  strain  lo  continue  propa¬ 
gating  along  the  filament  towards  the  impact  point  and  evertually  reflecting 
back  again,  and  a  certain  percentage  to  be  reflected  back  towards  the  fixed 
boundary  as  a  tension  pulse. 

Figure  6  illustrates  the  transverse  filament  displacement  history  asso¬ 
ciated  with  the  longitudinal  strains  described  in  Figure  5. 

The  above  computations,  even  though  based  on  an  approximate  linear 
stress-strain  behavior,  show  remarkable  agreement  with  the  experimental 
results  (such  as  those  shown  in  Figure  1,  for  high  tenacity  nylon  yarn). 
There  is  agreement  on  the  shapes  of  the  transverse  displacement  wave  be¬ 
tween  Figure  1  and  Figure  6.  TVds  is  evident  in  the  configurations  of  the 
filament  after  successive  reflections  of  the  transverse  wave,  as  indicated 
m  the  261  and  391  psec  displays  and  the  bend  in  the  130  psec  curve  at 
X=2.5  cm,  which  is  caused  by  substantial  differences  in  tension  level  across 
the  longitudinal  wave  front  at  that  position  (see  Figure  5).  In  addition 
there  is  quantitative  agreement  in  the  measured  and  predicted  velocities  of 
the  leading  edge  of  the  transverse  displacement  wave  and  with  the  projectile 
energy  losses  discussed  below. 

Explanation  of  the  Experimental  Results  for  Projectile  Energy  Loss 

An  explanation  of  the  experimental  curve  of  projectile  energy  loss  as  a 
function  of  striking  velocity  in  Figure  2  can  be  obtained  by  applying  the 
above  results  to  resolve  the  amount  of  kinetic  and  strain  energy  stored  in 
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<f  Of  IMPACT 


DISTANCE  cROM  IMPACT  X (cm) 

Figure  5.  LONGITUDINAL  STRAIN  WAVE  FOR  VARIOUS  TIMES  AFTER  IMPACT  OF  FILAMENT. 
Impact  velocity  =  334  meters/sec.  Materiel:  Linear  approximation  of  high  tenacity  nylon  of  Ref.  4. 

E  (modulus)  =*  80  grams/denier. 


£df  IMPACT 


DISTANCE  FROM  IMPACT  X(cm) 

Figure  6.  TRANSVERSE  DISPLACEMENT  OF  FLEXIBLE  FILAMENT  VERSUS  DISTANCE  FROM  IMPACT 
POINT  FOR  VARIOUS  TIMES  AFTER  IMPACT.  Impact  velocity  «  334  meters/aec,  Material:  Linear  approximation 
of  high  tenacity  nylon  fiber  of  Ref.  4,  E  (modulus)  »  80  grams/denier. 
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the  fiber  system.  The  shape  of  the  curve  is  the  result  of  either  the  strain 
energy  or  the  kinetic  energy  of  the  fiber  playing  the  dominant  role  at  dif¬ 
ferent  impact  velocities. 

It  is  convenient  to  consider  three  velocity  regimes  over  which  the  means 
of  obtaining  material  properties  and  the  manner  in  which  the  fiber-system 
responds  (i.e.,  the  importance  of  wave  propagation  effects  and  boundaries) 
are  quite  different: 

Regime  I  corresponds  to  a  quasi-static  fiber  deformation  where  multiple 
wave  reflections  have  averaged  out  the  velocity  and  strain  fields  and  wave 
propagation  effects  are  negligible.  The  material  properties  for  these  rates 
of  deformation  are  most  conveniently  obtained  by  techniques  outlined  in 
Reference  3.  The  transverse  velocity  distributions  are  fairly  linear  from 
zero  at  the  fixed  boundary  to  the  value  of  the  projectile  velocity  at  the 
inpact  point.  At  low  velocities  the  kinetic  energy  of  the  fiber  is  negli¬ 
gible  compared  to  the  strain  energy  stored  in  the  fiber  at  break.  As  the 
velocity  approaches  zero,  the  curve  asymptotically  approaches  some  breaking 
strain  energy  associated  with  the  moderately  high  rates  of  loading  experi¬ 
enced  by  the  fibers  in  this  regime.  As  the  impact  velocity  increases, 
kinetic  energy  increases  as  the  square  of  the  velocity  and  becomes  a  major 
part  of  the  energy  stored  in  the  fiber  system. 

Regime  II  corresponds  to  a  truly  dynamic  system  response  in  the  sense 
that  tne  wave-boundary  interactions  cause  the  peaking  energy  absorption 
effect  shown  in  Figure  2.  The  material  properties  for  these  rates  of  de¬ 
formation  are  most  conveniently  obtained  by  techniques  outlined  in  Reference 
4.  A  given  impact  velocity  causes  a  specific  amount  of  initial  strain  in  a 
given  fiber.  The  strain  builds  15),  upon  reflection  (as  illustrated  in 
Figure  5),  until  the  breaking  strain  is  reached  after  a  certain  period  of 
time.  If  this  time  corresponds  to  the  time  it  takes  for  the  transverse 
deflection  wave  to  just  reach  the  boundary  (e.g.,  between  217  and  261  nsec 
in  the  situation  illustrated  in  Figure  6),  then  the  condition  of  maximum 
fiber  kinetic  energy  (the  entire  fiber  having  a  transverse  velocity  close 
to  the  impact  velocity)  and  the  condition  of  maximum  strain  energy  coincide, 
giving  rise  to  the  peak  energy  absorption.  At  higher  impact  velocities  the 
breaking  strain  is  attained  before  the  entire  fiber  is  involved  in  trans¬ 
verse  motion,  so  again  the  kinetic  energy  contribution  drops  off. 

kegime  III  corresponds  to  a  purely  dynamic  fiber  response,  in  the  sense 
that  tne  initial  strains  are  close  to  the  breaking  strains  and  the  fiber 
breaks  almost  upon  impact  before  any  wave  disturbance  reaches  the  bound¬ 
aries.  Rate-dependent  effects  play  a  major  role  in  the  deformation  mechan¬ 
isms  at  these  high  rates.  Considerations  in  Reference  5  and  the  appendix 
must  be  employed  to  arrive  at  a  reasonable  representation  of  the  rate-de¬ 
pendent  material  behavior.  Since  less  and  less  fiber  is  involved  in  de¬ 
formation  the  energy  absorbed  decreases  to  zero. 
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Impact  of  Filament  with  Non-Linear  Material  Properties 

Figure  7  illustrates  typical  results  for  the  longitudinal  strain  distribu¬ 
tion  in  a  fiber  with  a  non-linear  stress-strain  curve  subjected  to  transverse 
impact.  Comparison  with  analytical  results  of  Smith,  et  al.,1**5  are  also 
shown. 

The  initial  steep  wave  front  of  the  1%  strain  is  due  to  the  fairly  linear 
nature  of  the  T-e  curve  up  to  1%.  Note  that  the  finite  element  solution 
(open  circles)  leads  the  analytical  predictions  because  of  a  slightly  higher 
value  of  initial  modulus  used  in  the  finite  element  analysis.  The  curved 
portion  is  due  to  the  gradual  decrease  of  modulus  with  strain  amplitude  and 
the  steep  rise  to  the  plateau  is  due  to  the  linear  nature  of  the  T-e  curve 
at  higher  strain  levels.  There  is  also  a  slight  dip  in  the  T-e  curve  which 
might  mean  that  a  weak  shock  wave  has  been  formed  in  this  steep  rise  to  the 
plateau  region. 

Impact  of  Filament  on  Elastic  Foundation 

Figure  8  illustrates  the  longitudinal  strain  caused  by  transversely 
impacting  a  linear  elastic  fiber  system  backed  up  by  an  elastic  foundation. 

The  particular  foundation  stiffness  chosen  for  illustration  is  quite  high 
(K=50,000  gr^s/cm  of  fiber  length/cm  of  deflection)  and  represents  an  ex¬ 
treme  siturrir'n.  The  fiber  properties  are  again  the  linear  approximation 
of  high  tenacity  nylon.  Comparison  should  be  made  between  Figure  8  and  the 
results  in  Figure  5  for  the  same  fiber,  without  backup,  impacted  at  the  same 
velocity.  Experiments1  have  shown  that  the  breaking  strain  for  this  fiber 
is  approximately  12.5%.  The  strain  concentrations  (Figure  8)  near  the  impact 
point,  resulting  from  the  restraining  action  of  the  foundation,  cause  the 
fiber  to  reach  breaking  strains  after  only  38  usee.  The  strain  history  in 
the  fiber  without  backup  in  Figure  5  extends  to  only  130  usee  after  impact. 
Analysis  of  the  subsequent  strain  history  predicts  that  the  breaking  strain 
of  12.5%  is  attained  after  230  usee  and  that  the  strain  distribution  along 
the  fiber  at  that  time  is  fairly  uniform:  11%±1.5%.  The  average  strain  at 
break  in  the  fiber  is  less  than  4%  (see  Figure  8) .  Accordingly,  the  strain 
energy  stored  in  the  backed  fiber  is  less  than  15%  of  that  stored  in  the 
fiber  without  a  foundation.  The  kinetic  energy  of  the  backed  fiber  is  also 
less  than  15%  of  the  unrestrained  fiber  kinetic  energy  at  break  since  less 
than  2  cm  of  fiber  is  involved  in  transverse  aotion  at  38  usee  and  the  entire 
fiber  is  involved  at  230  usee  (see  Figure  6). 

This  marked  change  in  performance,  when  the  deflection  of  the  textile 
fiber  is  restricted,  indicates  that  there  is  a  need  to  study  the  consequences 
of  the  constraining  effect  of  the  body  on  the  in-place  armor  system  and  that 
there  might  be  a  need  to  specify  some  deflection  criteria  for  flexible  armor 
systems . 

This  casv;  also  illustrates  the  versatility  of  the  analysis  to  handle 
easily  resistive  forces  which  cause  horrendous  mathematical  complexities  in 
an  analytical  solution. 
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Figure  7.  LONGITUDINAL  STRAIN  DISTRIBUTION 
IN  A  FIBER  HAVING  NON-LINEAR  STRESS- 
STRAIN  CURVE,  a.  Measured  stress-strain  curve  for 
high  tenacity  polyester?  b.  Longitudinal  strain 
distribution  at  24  /tsec  after  imcact  and  at  an  impact 
velocity  of  308  M/sec.  Finite  dement  analysis  o  o  ». 
analytical  result —  . 


DISTANCE  FROM  IMPACT  Xfcnt) 

Figure  8.  LONGITUDINAL  STRAIN  DISTRIBUTION 
PRODUCED  BY  THE  IMPACT  IN  A  FILAMENT 
SUPPORTED  BY  AN  ELASTIC  FOUNDATION. 
Material  is  the  linear  app'oximstion  of  high  tenacity 
nylon4  with  E  =  80  gm/denier.  Impact 
velocity  =  334  M/sec.  Elastic  foundation  K 
(stiffness)  =  50  X  103  gm/cm? 


Transverse  Impact  of  a  Flexible  Membrane 

Figure  9  illustrates  the  radial  .aid  circumferential  stresses  produced  in 
a  membrane  or  2-dimensional  sheet  of  flexible  material  for  various  times 
after  transverse  impact. 

Note  that  circumferential  buckling  occurs  during  the  propagation  of  the 
stress  waves  prior  to  reflection  at  the  boundaries.  The  buckling  is  caused 
by  the  fact  that  material  is  being  fed  into  a  smaller  annular  space  at  a 
rate  which  is  faster  than  the  buildup  of  radial  tensile  stress  and  the  asso¬ 
ciated  circumferential  Poisson  contraction  which  is  required  to  take  up  this 
additional  material. 

The  subsequent  stress  history,  not  illustrated  here,  predicts  the  effect 
of  wave  reflections  from  boundaries.  The  boundaries  have  a  much  less  dra¬ 
matic  effect  than  they  do  in  the  fiber  problem  since  the  stress  and  strain 
amplitudes  have  decayed  considerably  with  distance  propagated.  However,  the 
reflected  stress  wave  does  bring  the  circumferential  stress  into  tension. 

This  is  the  first  ligorous  model  to  be  presented  which  can  predict  the 
transient  deformation  of  a  membrane  subjected  to  transverse  projectile 
impact . 
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Figure  9.  CALCULATED  STRAIN  DISTRIBUTION 
IN  AN  IMPACTED  FLEXIBLE  MEMBRANE  USED 
AS  AN  APPROXIMATION  TO  THE  RESPONSE 
OF  A  NYLON  TEXTILE  FABRIC.  Nylon  material 
properties:  Young's  modulus  E  =  40  X  106  gm/cm? 
Poisson's  ratio  V  -  0.35;  density  p  =  1.15  gm/cm? 
thickness  7}  =  0.005  cm.  Impact  velocity  =  100  M/sec. 

a.  Radial  stresses  versus  radial  distance  from  impact. 

b.  Circumferential  stress  versus  radial  distance  from 
impact. 


CONCLUSIONS 

A  "finite  element"  mechanical  model,  capable  of  utilizing  experimentally 
obtained  material  prop  rties,  has  been  developed  for  analyzing  the  behavior 
of  a  constrained  fiber  system  when  subjected  to  transverse  impact. 

The  model  has  been  applied  successfully  to  the  following  transverse  im¬ 
pact  problems:  (1)  the  prediction  and  interpretation  of  experimental  re¬ 
sults  of  transversely  impacted  finite  lengths  of  nylon  yams  rigidly  clamped 
at  both  ends;  (2)  the  impact  of  infinitely  long  fibers  with  non-}*' near  ma¬ 
terial  properties  and  the  comparison  of  results  with  published  theory;  (3) 
predictions  of  strain  concentrations  caused  by  an  elastic  foundation  trans¬ 
versely  supporting  a  fiber  subjected  to  impact;  and  (4)  the  transient  de¬ 
formation  of  an  isotropic  membrane  subjected  to  projectile  impact. 

In  view  of  the  versatility  of  the  finite  element  analysis,  as  indicated 
in  the  filament  applications  cited  above,  it  appears  reasonable  to  expect 
that  further  effort  will  yield  improvements  in  modeling  the  physics  of  an 
actual  fabric  system  and  ultimately  provide  a  viable  model  for  predicting 
the  response  of  the  in-place  anisotropic  fabric  or  felt  armor  system  sub¬ 
jected  to  transverse  impact. 
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Appendix.  LINEAR  VISCOELASTIC  EFFECTS 


In  the  analyses  outlined  above,  the  material  , j^perties  were  assumed  to 
be  rate  independent.  In  one  of  the  latest  papers  by  Smith5  the  strain 
distribution  along  the  fiber  was  measured  experimentally.  The  marked  dif¬ 
ferences  illustrated  in  Figure  A-l  between  these  results  and  the  analytical 
results  such  as  depicted  in  Figure  7  of  the  main  text  were  attributed  (ac¬ 
cording  to  an  excellent  physical  argument)  to  viscoelastic  effects. 

The  finite  element  model  was,  therefore,  modified  to  accept  a  general 
type  of  stress  relaxation  modulus  which  can  be  expressed  analytically  as: 

Ect-o  A1 

This  convolution  integral  was  integrated  by  finite  differences  in  order 

to  be  cast  into  a  form  amenable  to  the  finite  element  model  analysis  follow¬ 
ing  the  procedure  in  Reference  8.  Integrating  the  equation  by  finite  dif¬ 
ferences  in  equal  time  increments  t(i)  (i=l,2,3, . . .  ,r.+l)  where  t(l)— 0  and 
t(n+i)=t  and  using  the  notation 

T(n+1)  =  T  (A,t(n+1)) 

£(J,n+l)  =  E(x,t(n+1))  A2 

E(n+i)  =  E(t(n+1)) 
the  stress-strain  relations  become: 

T(0=£{i>ai) 

4  A3 

Tu)  =  e(2 ){ i  [  ECO  *  ec  2)] }  -  euKHsco-  > 

The  only  modification  required  injthe  analysis  is  to  keep  track  of  the 
strain  history  of  each  element  E(x,n+1)  and  compute  the  current  stress  ac¬ 
cording  to  the  above  equations  (A3) . 

The  problem  of  transversely  impacting  the  linearly  elastic  filament  as 
shown  in  Figure  5  was  reconsidered,  using  a  linear  relaxation  modulus  vary¬ 
ing  from  100  gm/ denier  to  60  gm/ denier  in  40  usee  as  shown,  in  the  insert  to 

Figure  A-2. 

Figure  A-2  shows  the  results  of  the  rate- independent  finite  element 
model  vs  the  linear  viscoelastic  model.  The  particular  relaxation  time 
used  was  based  on  some  of  the  qualitative  arguments  by  Smith  et  al.s 
The  results  qualitatively  agree  with  the  results  of  Smith  et  al  illustrated 
in  Figure  A-l  and  tend  to  prove  his  contention  that  there  is  a  distinct  re¬ 
laxation  process  on  the  order  of  50  usee  in  the  subject  nylon  fiber. 

A  slight  modification  of  the  program  to  account  for  non-linear  stress 
strain  curves  exhibiting  relaxation  behavior  is  now  needed.  These  modifica¬ 
tions,  along  with  more  judicious  interpretation  of  the  ’'rate-independent" 
moduli  obtained  by  Smith  et.  al.,4  should  provide  a  model  capable  of  pre¬ 
dicting  very  accurately  the  detailed  strain  distribution  (Figure  A-l)  in 
textile  fibers  caused  by  transverse  impact. 
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a.  Stress-strain  curve  for  high  tenacity  nylon. 


b.  Rate-independent  prediction  obtained 
from  stress-strain  curve  (solid  line). 
Experimentally  measured  strain 
distribution  (dashed  line) 


Figure  A  1.  COMPARISON  OF  LONGITUDINAL  STRAIN  IN  FILAMENT  SUBJECTED  TO  TRANSVERSE 
IMPACT  CALCULATED  FROM  RATE-INDEPENDENT  THEORY  AND  MEASURED  EXPERIMENTALLY? 
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a.  Linear  relaxation  model  representing  a 
viscoelastic  approximation  of  high 
tenacity  nylon? 


b.  Strain  distribution  obtained  from 
rate-independent  theory  and  from 
the  viscoelastic  approximation 
(curve  represents  the  experimental 
behavior  depicted  in  A-lb). 


Figure  A-2.  LONGITUDINAL  STRAIN  IN  VISCOELASTIC  FILAMENT 
SUBJECTED  TO  TRANSVERSE  IMPACT 
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A  "finite  element"  mechanical  model,  capable  of  utilizing  experimentally  deter¬ 
mined  material  properties,  is  presented  for  analyzing  the  deformation  processes  in  a 
finite,  constrained  flexible  filament  subjected  to  transverse  impact. 

The  model  is  then  used  to  analyze  the  results  of  experimentally  (transversely) 
impacted  finite  lengths  of  nylon  yams  rigidly  clamped  at  both  ends.  The  experimental 
analysis  is  presented  in  Part  I  of  this  work. 

The  analysis  is  used  to  determine  the  effects  of  an  elastic  foundation  trans¬ 
versely  supporting  a  flexible  filament  subjected  to  transverse  impact. 

Finally,  results  of  a  finite  element  analysis  of  a  transversely  impacted  iso¬ 
tropic  membrane  are  presented  to  indicate  the  direction  for  future  research. 

An  appendix  is  attached  which  demonstrates  the  ability  to  include  rate -dependent 
material  properties  in  the  analyses  described  above.  (Author) 
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